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1 INTRODUCTION 

In this talk we review some recent results from random matrix models as 
applied to some non-perturbative issues in QCD. All of the issues we will 
discuss touched upon the important phenomenon related to the spontaneous 
breaking of chiral symmetry. 

Instead of using the standard lore of Green's functions in random matrix 
models, we will insteacLchoose to work with their functional inverse or Blue's 
functions as defined byEJ. This way of doing things sheds much insight into the 
physics of random matrix models, and is probably the most "natural" and user- 
friendly way of introducing the formal, but powerful mathematical concept of 
free random variables In this talk we skip almost all mathematical details 
(referring the interested reader to the published work) , and append some formal 
and mathematically relevant issues to an Appendix. 

Exploring new and pertinent mathematical issues in the realm of frontier 
physics has been one of Mannque Rho specialty. We hope that the present 
lecture which combines topics of current physical interest in chiral symmetry 
breaking and QCD, with novel mathematical concepts, can be regarded as a 
small tribute to Mannque Rho's 60th birthday. 

2 Chiral Symmetry Breakdown and Random Matrix Models 



"Talk presented by MAN at "Hadrons in Dense Matter", Seoul, 1997. 



1 



2.1 Disorder 



The first important question is: why random matrix models have anything in 
common with QCD? The answer will come from the generic aspects of the 
way chiral symmetry is spontaneously broken in the vacuum, and the fact 
that sigma-models are good starting point. In the process, we will show that 
the results of RMM fall into universal (microscopic regime) and non-universal 
(macroscopic regime) both of which give us much insights to the intricacies 
of the QCD Dirac spectrum in a finite volume. We would like to stress from 
the onset that the universal results provided by RMM are not only applicable 
to QCD but to any effective model of QCD that breaks spontaneously chiral 
symmetry. 

In the macroscopic regime qualitative but hopefully generic features of 
non-perturbative results can be probed with much insights to the important 
parameters at work in reaching the thermodynamical limit. QCD in a finite 
volume as probed by current lattice simulations is in many ways a complex 
and disordered system. RMM are attractive tools to model disorder in the 
presence of external and deterministic parameters such as masses, temperature, 
chemical potential, and vacuum angle. Changes in these parameters introduce 
specific signatures on the QCD Dirac spectrum, some of which may be generic 
enough to be modeled by pertinent matrices. As most of the physics retained 
in RMM is soft, the crucial assumption relies on the decoupling of soft and 
hard modes in QCD. The many qualitative agreements between RMM results 
and lattice simulations suggest that this may be the case. A thorough analysis 
can be achieved by thorough lattice cross-checks between cooled and uncooled 
simulations. 

Now, consider the Banks-Casher relatione! 

This remarkable relation links the order parameter of chiral symmetry break- 
down, i.e. quark condensate (qq), to the average spectral density v of the 
eigenvalues near A = of the massless Dirac operator ilp(A). The presence 
of the Euclidean volume V4 in (|l]) is not accidental. It points out, that the 
average spectral density of the Dirac operator at zero virtuality has to grow 
like volume V4, in order to provide the mechanism for non-zero condensate in 
the thermodynamical limit. Free quarks in a finite volume V exhibit at zero 
virtuality a density \fV^ which is not enough in the thermodynamical limit to 
cause the spontaneous breaking of chiral symmetry. 

We note the similarity of Banks-Casher relation to Kubo-Greenwood for- 
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mula for direct current conductivity a 

o = Dp(Ep) (2) 

relating the conductivity to density of electronic states at the Fermi level, via 
the diffusion constant D. This analogy suggests, that the regime of sponta- 
neous breakdown of chiral symmetry is actually diffusive and originates from 
delocalization of quark eigenmodes. A thorough characterization of this regime 
and others may be found in our recent worka. 

In real QCD, the soft modes are likely to set in around instantons config- 
urations or any of their topological relatives, e.g. monopoles. Cooling of the 
lattice has been shown to reproduce the gross features of QCD in the infrared 
limit, so a simplification of the quark modes to the ones near the zero virtuality 
regime may be qualitatively legitimate for bulk observables and Dirac spectra. 
This is minimally achieved through the use of RMM. The latters capture the 
essentials of the quark zero modes in a random instanton vacuum El, by which 
they were originally inspired. 



2.2 Microscopic universality in vacuum 

The fact that the spontaneous breakdown of chiral symmetry along with its 
explicit breaking provides powerful constraints on on-shell amplitudes is widely 
known and forms the basis of chiral perturbation theory!. What is less known 
is that the same mechanism provides a severe constraints in the microscopic 
regime. To show this, consider the pion propagator in a periodic Euclidean 
box, 



1 ^-^ e ~iK n x ^ 1 e ~iK n x 
Gtt(x) = — 2_] "772 i 2 = 77 + T?2l^2 _ 2 ( 3 ) 



where the primed sum is over non-zero modes 2irn,i/L, (i = 1,2,3,4). In the 
macroscopic regime mV ^> 1, Gell-Mann-Oakes-Renner relation m\ = m£/F^ 
(where S = allows the usual chiral power counting, leading to standard 

chiral perturbation theory. 

In the microscopic regime mV <C 1, the primed sum is negligible, since 
the zero mode contribution (mV)^ 1 dominates (||). The chiral counting breaks 
down. However, it is still possible taxesum the zero modes. The general result, 
obtained by Gasser and LeutwylerQ could be rewritten in terms of partition 
function: 

Z(m) = [ dUe v ^ mU+h - c - ] (4) 
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Note that partition function does not involve any derivatives. The partition 
function depends solely on: 

(i) - The nature of the coset (pattern of spontaneous breakdown of chiral sym- 
metry), hidden in the integration measure. 

(ii) - Explicit breaking of chiral symmetry ( representation (Nf,Nf).) 

Next important observation, due to Leutwyler and Smilgai, is based on 
the following transformation based on the 9 angle 

Z n {m) = 1^ — e - m9 Z(me ie / N ' ) = / dU{det U) n e v( - mU+h -^ (5) 

JO 27T J 

L.h.s., by definition, describes the QCD partition function for fixed topological 
(winding) number n. R.h.s. is rephrased in hadronic variables. Comparing 
the same powers of m on both sides yields an infinite hierarchy of sum rules, 
relating the eigenvalues coming from quark determinant in Z n (m) to hadronic 
parameters on the r.h.s. of (|5|). For example, comparison of quadratic terms 
leads to 



vA^Xl/ n A(n + N f ) 



(6) 



k 

It was conjectured by Shuryak and Verbaarschot i, that there must exist a 
master formula, based on some microscopic spectral density v s (x), where x = 
VX, which generates, via taking suitable moments, the LS sum rules. 

At this step the impatient reader may ask, how do the random matrices 
come into this scenario? The answer is simple - let us construct the simplest 
(minimal) model which satisfies both constraints (i) and (ii). For pedagogical 
reason let us constrain temporarily to one flavor. Such a model is a four-fermi 
model in dimensions: 

L = qUnuP* 9 q + yq R q L q R q L (7) 

where the dimension of Grassmann fields is ul, ( n R) f° r "left" ("right") com- 
ponents, respectively, analog of 75 is defined as (l nR ,— 1„ L ) and the "volume 
V" of the system is nothing but til + Ur. It is not surprising at all, that such 
a model is consistent with restrictions of chiral symmetry breaking, since, it is 
just the zero mode sector of usual Nambu-Jona-Lasinio model! 

What is perhaps more surprising, is that this model can be_r.ewritten in 
terms of the bosonic variables A ab = qlq 3 ^ in the following forma 

ime' 



Z(9,m) - / dAeT VVM ' det ( "J t ,J o _ ie ) (8) 
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But this is, by definition, the partition function for the chiral (cf. block struc- 
ture), random Gaussian (cf. quadratic potential in the exponent) matrix (cf. 
matrix-valued field A) model. 

This identification of the class of random matrix models allows to use 
the whole arsenal of RMM to get the microscopic spectral density—Using the 
variant of orthogonal polynomial method, Verbaarschot and ZahcdEH found an 
explicit form of microscopic spectral density for massless QCD 

u s {x) = ±Y?x{J% f {Vx) - J 2 Nf+1 {Vx)J 2 Nj _ x {Y.x)) (9) 

It is a straightforward exercise to show, that the moments of microscopic 
spectral density generate, via doable integrals over Bessel functions, the r.h.s. 
of the Leutwyler-Smilga sum rules. 

The last unresolved issue is the universality, i.e. the independence of the 
microscopic spectral density on the form of the matrix potential (one could 
easily visualize higher powers of AA' in the exponent). At first sight, such 
universality looks like a miracle, and only very recently was prcsien to hold for 
arbitrary series-like potentials by Damgaard and collaborators using rather 
formal methods of RMM. Actually, the physical reason of universality is simple 
- all higher order than quartic terms in NJL-like lagrangian are suppressed in 
the thermodynamical limit (i.e. behave like 1/yV). Let us finally mention, 
that sum rules and microscopic spectral densities were successfully checked in 
several models of the QCD vacuum, e.g. instantpn liquid model E3. They were 
also checked with high accuracy on the lattice Ej, albeit in the last case the 
form of microscopic universality is of different type (for SU(3) C Kogut-Susskind 
fermions), due to the known ambiguities caused by Nielsen-Ninomiya theorem 
(fermion doubling). 

2.3 Microscopic universality in matter 

We may ask, to what extent the microscopic universality survives the effects 
of matter, e.g. the temperature or presence of the baryonic potential. Generic 
arguments from previous section based on power counting suggest that most 
of the observations survive the changes in external, deterministic parameters, 
provided that chiral symmetry is still broken. To make this argument more 
transparent, let us bring the standard universality argument by Pisarski and 
Wilczek E£l based on mapping of the QCD onto the linear sigma model. For 
exemplary case of two-flavor QCD in the vicinity of phase transition, the a 
model takes the form 

L = Tr(m$ + h.c.) + g (T)Ti&$ H (10) 
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with go(T) proportional to the reduced temperature. The only effect for the 
crucial for LS-sum rules first term in this Lagrangian is to "dress" thermally 
constant modes: $ — > <£>(T). We therefore suspect that all LS sum rules 
and microscopic spectral densities will stay unchanged, provided we replace 
the scale £ of RMM by S(T). The explicit form of temperature dependence 
of the scale is not universal, i.e. depends on the potential of RMM. Here, 
equivalently, depends on the explicit form of relevant operators (dotted terms) . 
This generic scenario agrees with detailed calculations Ej. Similar arguments, 
modulo subtleties related to the formation of the nucleon Fermi surface, could 
be presented in the case of baryonic potential. 

Does it mean that the micposcopic spectral formula is always sacred? Not 
necessarily. It was pointed inllJ that when one also rescales the masses, i.e. 
in the limit y = mV ~ Ay = x, microscopic spectral density will become 
function of two variables, i.e. vs{x,y), provided unqucnching of the fcrmionic 
determinant. The universality of-these "double microscopic" spectral densities 
was proven by Damgaard et al. 113 Despite the double limit procedure is not 
physical, the new sum rules generated by double spectral densities could serve 
as a theoretical concept allowing better assessment of spontaneous breakdown 
of chiral symmetry on the lattice. 



2.4 Structural changes of Dirac spectrum - temperature 

In this subsection we will try to get some insight how external parameters may 
influence the spectrum of random matrix model. The influence on microscopic 
features was discussed above. Here we would like to model the non-universal 
(c/. $(T)), albeit perhaps qualitatively correct thermal behavior of quark con- 
densate. Let us start again from Banks-Casher relation. For zero temperature 
it is easy to check, that indeed the average spectral density of RMM is non 
zero for A = 0. We expect it intuitively anyhow, since we know that RMM is 
just the longest wavelength limit of NJL model, but let us recover this result 
using RMM point of view. 

In practice, it is more convenient to work with Green's functions then 
spectral densities. Let us therefore define the Green's function for Gaussian 
random matrix model (for 9 = 0): 

For our purposes (zero virtuality) we need to calculate this resolvent at z = 
only, or equivalently, for any z with m = and finally put z = im. Here 
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matrices A are hrXul. The two diagonal blocks have implicit size Uri^ xn^j. 
For fixed n R / Ll we use the chiral basis 

G = <7l + <fc75 (12) 

with 1 = diag(l ni5 ,, l nL ) therefore, after performing the averaging (..) with 
Gaussian measure, we get chirality even and chirality odd resolvents 

TrG = g + xg 5 = ~(z - ^ z 2 - 4 + 4x 2 / z 2 )\ z = lm 

X 

TrG7 5 =xg + g 5 = —\ z =%m (13) 

Since spectral densities are just related to imaginary part of the Green's func- 
tions^], we immediately read out relevant spectral densities 

u + (X, x) = \x\6(X) + -L^^-A^-A*) 

v-(X,x) =xS(X) (14) 

with Aj_ = 2 ± 2yi — x 2 and asymmetry x = (tir — ni)/2N . Note that for 
x = one recovers Wigner's semicircle lawEi In particular, the chirality even 
spectral function v+ does not vanish at A = 0, leading therefore, via analogy 
to Banks-Casher relation, to "quark condensate" . 

The chirality odd part V- is a direct measurement of the difference in 
the spectral distribution between left-handed and right-handed "quark zero 
modes" . If resembles naively Atiyah-Singer theorem, despite here the origin of 
the asymmetry is purely kinematical ( "rectangularity" of the ensemble), and 
not related to any topology. We will come back to this point when discussing 
the C/(l) anomaly, so in the mean-time let us put x = 0. 

Now we see how to introduce external parameters into the random matrix 
model - all we have to do is to bosonize NJL model in matter and look at the 
softest modes. Therefore the modifications of the medium correspond to to 
the modification of Q by 

/ ifi-/z\ / m iA\ , , 

\iSl + fi ) + [iAl m ) 

where \i is the chemical potential for "quarks" and fi = u> n l n ® Ijv- Here 
uj n are all fermionic Matsubara frequencies. For simplicity, we restrict our- 
selves to the lowest pair of Matsubara frequencies ±7rT and the chiral limit. 

6 Note l/(y + £0) = PVl/y - mS(y). 
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This corresponds to the model proposed inEl To investigate chiral symmetry 
breaking in this model we should calculate the Green's function and, through 
Banks-Casher relation, obtain the condensate. We have therefore to "add" , 
modulo trivial overall factor i, schematic deterministic (D) and random (R) 
"Hamiltonians" . 



nT\ f A 
ttT I + \ At 



(16) 



At this moment we could use some of the Blue's magic, here addition theorem. 
The Blue's function for the model (|l6| ) satisfies (see Appendix) 

B(z)-z = B D (z) (17) 

where Bo is the Blue's function of the deterministic piece. By definition we 
have 

z - G(z) = B D (G(z)) (18) 

Now if we evaluate the Green's function of the deterministic piece on both 
sides of this equality we will get so-called Pastur-Wegner equationcil 

G D {z - G(z)) = G D {B D {G(z))) = G(z) (19) 

The explicit form of deterministic Green's function is Gd{z) — ^( z _ 1 7rT + 
z+Vt ) ' smce the deterministic eigenvalues are ±7rT (lowest Matsubara fre- 
quencies), so the final result reads 

G 3 - 2zG 2 + {z 1 - ir 2 T 2 + 1)G - z = (20) 



Figure 1 shows the imaginary part of the relevant solution of (|20|) as a function 
of eigenvalues and temperature. || The structural change at the forking point 
corresponds to vanishing i>{\ — 0), therefore restoration of chiral symmetry. 



2.5 Critical temperature 

Critical temperature, corresponding to the point of restoration of the "chiral 
symmetry" of this model, could be inferred from differentiation of the Blue's 
function corresponding to (^o|). Intuitively it is obvious - the endpoints are 
branching points, so the derivative of the Green's functions blows up, hence 
the functional inverse of Green's function (Blue's function) takes extremumEI. 

c The other two solutions are either real or lead to negative spectral density. 
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The condition B'(G) = 0, together with ( |2C| ) leads to the real endpoints of the 
spectra located at 



1 (47r 2 T 2 -l + V87r 2 T 2 + l)i 
V&nT V8tt 2 T 2 + 1 ~ 

1 (4tt 2 T 2 - 1 - V8tt 2 T 2 + 1)1 
V8ttT V8tt 2 T 2 + 1 + 1 

The condition I2 = 0, corresponding to the situation when two-arc support 
of the spectrum [— Ii,— 12] U [/2,-fi] evolves into the single interval on the cut 
[— 7i,+Ji] corresponds to "phase transition" and defines the critical tempera- 
ture T» = I/71", in units where the width of the random Gaussian distribution 
is set to 1. 



±h = 
±h = 



2.6 Critical exponents 

The set of all critical exponents could be easily inferred from solving algebraic 
equation ( ptj|) and calculating the free energy of the system through integrating 
Blue's functions (see Appendix) giving free energy 

^ = G 2 +logi__^ (22) 

Explicitly, they are: (a, 0, 7, 6, v, rj) = (0, 1/2, 1, 3, 1/2, 0), respectively. There 
are mean-field type independently on the piimber of Matsubara frequencies 
used, in agreement with general arguments E3. 



2.7 Chemical potential - "phony vacua" 

Let us now move to the case of finite chemical potential, switching the tem- 
perature off for the clarity of the presentation. The matrix model is now given 
by 

'0-/A (0 At' 
H J \A 



(23) 



and was originally solved by StephanovE3 using standard RMM tools. The 
deterministic part is now non-hermitian (eigenvalues are complex numbers) and 
we have two distinct regions (at least in the quenched case). The Green's 
function in the holomorphic region outside the blob of eigenvalues is given by 
@ with the substitution tt 2 T 2 — ► -/i 2 : 

G 3 - 2zG 2 + (z 2 + /i 2 + 1)G - z = (24) 
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Figure 2: Schematic mapping illustrating conformal transformation of the cuts supporting 
the eigenvalue distribution visible at figure 1 onto the boundaries of islands for complex 
eigenvalue distribution. The shaded regions represent holomorphic domain, where the map- 
ping is valid. 



Instead of following the original argumentE3, we suggest again the use of Blue's. 
In order to find the boundary of the nonholomorphic region we may exploit 
the method of conformal mapping via the Blue's functions (see schematic fig- 
ure 2), and transform the cuts of the T ^ case into the boundary by the 
transformation 

z — >w = z- 2G(z) (25) 

where G(z) is the appropriate branch of the cubic equation ( p0[ ) with formal 
identification tt 2 T 2 = +[i 2 . 

The result of this mapping (continuous line) is presented on Fig. 3. The 
dots are the numerically generated eigenvalues of (|2^) . The mapping, repro- 
duces immediately results of original random matrix model for fj, E2I. The 
islands reproduce the domains of mixed-quark condensate ("phony vacua"), 
obtained as an artifact of neglecting in the process of averaging the phase of 
the complex fermion determinant. This condensate is not the usual chiral quark 
condensate, but quark-conjugate quark condensate which forms due to the fact 
that quenching of the phase stops penalizing attraction leading to such con- 
figurations. Such a condensate violates spontaneously baryon number.^] This 
model suggests the reasons for failure of the quenched lattice calculations with 

d Vafa-Witten restrictions do not hold for complex measure. 



11 




-2-1012 -2-1012 -2-1012 

XXX 



Figure 3: Comparison of the numerically obtained distribution of eigenvalues for few sample 
chemical potentials with analytical results. The solid lines represent the boundaries obtained 
by conformal mapping. 



chemical potential, due to the appearance of this mixed-condensate or, equiva- 
lently, infinitely growing fluctuations when approaching the boundaries of the 
islands E3, due to the appearance of so-called "baryonic pion". 



2.8 Chemical potential - ordinary vacua 

We may also consider the random model without suppressing the phase of 
the determinant. The investigation of the behavior of uaquenched partition 
function (e.g. location and densities of Lee- Yang zeroesEJ) could be immedi- 
ately inferred from integration of Blue's functions. The location of singularities 
separating the "phases" of the system may be found from equating the real 
parts of the free energies calculated for different solutions (labeled of the 



cubic Pastur-Wegner equation 



partition function for model (|22 ) is given by 



24) . Explicitly, the location of singularities of 



Z G j . . . , Z — G j 



ReFi = Re(G'i + log —t^——) = Re(G J 2 + log 1 ) = ReF 3 (26) 
Gi Gj 



This procedure describes very nicely precise (up to 500 digit accuracy) numer- 
ical calculations ofc3. The_&esulting phase structure suggests first order phase 
transition, as discussed in £3. 



2.9 Phase diagram 

The combination of temperature and chemical potential can be easily investi- 
gated as discussed in subsection 2.4. The 1/N approximation is reminiscent 
of the mean-field treatment discussed originally in 113. The phase diagram 
(for chiral case) for condensate is plotted on Fig. 4. For chemical potential 
^ > Mer it = 1/611 ln(2 + a/3) the transition is first order (due to lack of vector 
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Figure 4: Condensate vs temperature and chemical potential for chiral random matrix model 
with infinitely many Matsubara frequencies included. 

couplings in Od NJL model) and disappears at /i = 0.5S. At zero chemical 
potential, the transition is second order (mean-field) with critical temperature 
T C rit = 1/4E. The different numerical coefficient comparing to subsection 2.4 
stems from the fact, that here we have taken into account all Matsubara modes. 
It is interesting to compare Fig. 4, obtained using only the constant modes, to 
the full NJL result of Weise and LutzEj and also Zhuang et al.EZI 

2.10 U(l) problem 

Let us now look at more subtle applications of random matrix models, where 
usual mean field arguments do not hold and the dynamical effects of fermionic 
determinant have to be taken into account. We start from the U(l) problem, 
i.e. the fact that pseudoscalar boson is much heavier than other pseudoscalars, 
violating bounds based on chiral symmetry. This problem is related to ax- 
ial anomaly. In subsection 3 we have encounter the "kinematical version" of 
Atiyah-Singer theorem. Since in QCD Atiyah-Singer theorem is nothing but 
the global (integrated) axial anomaly, we may hope that the effects of fluctua- 
tions in (jir — til), leading to chirality-odd spectral density, may mimic some 
features of anomalous sector of QCD. 
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Naively, it looks that the effects of asymmetry x are meaningless in the 
large N limit. Variable x — (tir — til)/2N seems to be irrelevant in the large 
TV limit. This reasoning is false, due to the fact, that the subleading at the first 
glance terms coming from asymmetry are enhanced by factor N coming from 
fermionic loops. The resolution of this problem requires therefore an infinite 
resummation of contributions originating from fermion determinant. 

Technically, we could do it in one step using the integration of the Blue's 
(see Appendix). For a fixed asymmetry x = x/2N we have 

d z log Z(x) = -2NN f iTrG (27) 

where TrG is given by first relation in (|l^). When integrating Blue's or Green's 
functions we need only to study the contribution of order x 2 from Green's func- 
tion, therefore G(z) w — x 2 / (z 2 V z 2 — 4). Straightforward integration leads to 

pogZ(x)] a = ±N f - x 2 (28) 

-^v rniymr + 4 + m) 

where index reflects the approximation and the constant was fixed from im- 
posing the right asymptotic behavior. Full partition function is 

Z = J dxe^Z(x) (29) 

where integration over \ takes into account the fluctuations in number of right- 
handed and left-handed zero modes. The Gaussian measure is consistent with 
axial anomaly E3. 

Inserting (E8h into above equation leads to replacement 



Nf _ N 



i/x* - i/x* = Vx* +E — r— — ) = Vx*+X> (so) 

i=i mf(Jm 4 f +4+m,f i=i 

where we have used the shorter notation and we have reinstated more general 
(unequal flavors) flavor dependence. 

Note that for Nf = (quenched) the topological susceptibility is just 
the given by x* j a $ it should, reminding the Witten-Veneziano relation known 
from QCD without light flavors. The unquenched result x* shows the screening 
caused by fermion determinant (quark loops) and vanishes in the chiral limit, in 
agreement with the theorems of the QCD. The formal trick with the integration 
of the Blue's functions allowed us to resum infinitely many infrared terms, and 
obtained a strictly non-perturbative result. "Quark loops" (terms of order x 2 ) 
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from fermion determinant were absolutely crucial for this resummation, since 
quenched result does not lead to the screening. The presented result is the 
analog of the screening phenomenon in a statistical ensemble made of negative 
and positive topological chargesJii full (unquenched) QCD. Similar mechanism 
works in the instanton vacuumed. 



2.11 Pseudoscalar susceptibility 

Let us finally demonstrate that the manifestations of the U(l) anomaly are 
identical in QCD and our random matrix model, leading to the same generic 
screening of pseudoscalar susceptibility. We note that the RMM leads to 
"Ward identity" relating topological susceptibility, pseudoscalar susceptibility 
and condensate 

2m m 2 
X = - W ^G- W X P s, (31) 

in analogy to anomalous Ward identity in Minkowski QCD 

un — ttx^ f — — 

iXtop = -jp(W>) + ^2 J d^xiT^i^x^i^iO)) (32) 

where the pseudoscalar susceptibility is defined by the integral in (|32| ), The 
resolution of the U(l) problem in QCD stems from the fact that for small 
chiral mass m, the absence of the U(l) Goldstone boson requires that \top = 
— m^ip)/^] to order 0{m 2 ), or in other words, pseudoscalar susceptibility 
has to be finite in the chiral limit. In random matrix model, the above un- 
quenched calculations and Ward identity yield the pseudoscalar susceptibility 
to be equall23 



Xps 



— [-(trQ 1 75)+"^(trQ X 75trQ 1 75) connecte 



(E-) 2 



4^7* (33) 



therefore finite in the chiral limit. The first equality is the definition of pseu- 
doscalar susceptibility in chiral random model. The second equality is the 
result of our calculations. Note the non-trivial cancellation of the infrared sin- 
gularities 1/m 2 between the two terms in (|3|). This remarkable cancellation is 
the result of competition between the trace and trace-trace terms in (|33|). Naive 
expectation based upon large N counting suggests wrongly that the trace-trace 
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term is subleading. Neglecting the trace-trace term leads to diverging pseu- 
doscalar susceptibility in the chiral limit, therefore to the "Goldstone pole" in 
U(l) channel, and consequently to the failure of the resolution of the anomaly 
problem. Like in real QCD, the unquenched treatment of large oscillations 
(©(./V 1 / 2 ) , where N commensurates with volume) in topological charge are 
needed for proper physics in the U(l) channel. We found rather remarkable 
that the simple matrix model without any further dynamical assumptions is 
able to illustrate the mechanism of screening consistent with our knowledge of 
the QCD. We do hope also that some of these results are relevant for lattice 
simulation for U(l) observables, taking into account that the resummation of 
infrared terms might be troublesome on the finite lattice. 



2.12 Finite 9 

QCD at finite vacuum angle 9 is not well understood. Despite experimental 
evidence suggest that 9 < 10~ 10 , the issue of potential CP breakdown of QCD 
poses a serious theoretical challenge. Since due to the U(l) anomaly, the 9 term 
may be traded from gauge fields to the quark mass matrix, some problems of 9 
term may be investigated using the chiral random matrix models, provided that 
the mass matrix diag (mi, ma, ...mN f ) is replaced by e l8 ^ Nf diag (mi, ma, ...mN f ). 

Saddle point analysis cil leads to the following unsubtracted free energy 
(for small masses) 



Nf 

F = J2(P"~ 1 °SP'~ 2 — (34) 



2=1 

x*v 

2N 



N f N, 



2=1 



— Pi 



where saddle point solutions of RMM are parameterized as pje l ^ e ^ Nf ~^ j \ Re- 
markably, this RMM result reproduces Witten equations E3, derived using dif- 
ferent arguments. For Nf — 3 the solution at 9 = ir is doubly degenerated, 
implying spontaneous breakdown of CP, provided the Dashen bound holds. For 
physical masses, however, spontaneous breakdown of CP seems to be excluded. 

Encouraged by these analytical results, we may probe numerically the un- 
quenched 9 problem in the framework of random_rnatrix models. At the left 



side of Fig. 2.12 we plot the numerical results Ed for free energy as a func- 
tion of vacuum angle. On the right side, we compare the similar plot fromLj, 
obtained for dynamical CP 3, model (plots for QCD do not exist due to the 
complexity of the problem). Both plots exhibit a kink at some critical value, 
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Figure 5: Free energy as a function of 9/tt for RMM (left) and CP 3 modelEa. Here A = N/V. 

suggesting first order "deconfining" phase transition. If this transition would 
persist for QCD, such picture may suggest a dynamical resolution of the 9 
problem. Namely, deconfmcment may always happen due to the fact that 
chromomagnetic monopoles would acquire electric charges 9/2ir, and would 
screen in consequence long-range color forces. Therefore the presence of con- 
fining phase forces 9 = 0. 

Such scenario requires that the kink visible on both figures migrates to 
zero in the thermodynamical limit. 

RMM allows another scenario consistent with the CP 3 data. High precision 
numerical studies in RMM Ell strongly suggest that the major factor driving 
the position of the kink seems to be rather the ratio N/V, where N refers 
to number of zero modes in a fixed volume V. Therefore the kink observed 
for dynamical models like CP 3 may be the result of subtle interplay between 
quark masses, finite size effects and accuracy in numerical assessments. 



3 CONCLUSIONS 

In this talk, we tried to give a quick overview of applications of RMM to var- 
ious domains of QCD, with the emphasis on the phenomenon of spontaneous 
breakdown of chiral symmetry. The main intention of this presentation was 
to exhibit the generic behavior based on fundamental physical premises. We 
therefore have skipped several important technical details and we have often 
referred to original papers. Meanwhile, we tried to stress the fact, that several 
"miraculous" at the first glance agreements of RMM and QCD, either at quan- 
titative level (microscopic regime) or at qualitative level (macroscopic regime), 
are based on firm and well understood principles of spontaneous breakdown of 
chiral symmetry. 
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We do hope that this way of looking at RMM in reference to QCD may 
be helpful at least for a non-expert in Random Matrix Models. 
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APPENDIX - FIVE FACETS OF BLUE'S 

The fundamental problem in random matrix theories is to fpd the distribution 
of eigenvalues in the large N (size of the matrix M.) limit LJ. One encounters 
two generic situations — in the case of hermitian matrices the eigenvalues lie 
on one or more intervals on the real axis, while for general non-hermitian ones, 
the eigenvalues occupy two-dimensional domains in the complex plane. In both 
cases the distribution of eigenvalues can be reconstructed from the knowledge 
of the Green's function: 

where averaging is done over the ensemble of N x TV random matrices generated 
with probability 

P (M) = - e - m " v{M ^ (36) 

For hermitian matrices the discontinuities of the Green's function coincide with 
the eigenvalue supports. In the simplest example of the Gaussian potential 
V{M) = M 2 for hermitian matrix A4, a standard calculation gives 

G{z) = \{z-Jz^i) (37) 

The reconstruction of the spectral function [i.e. eigenvalue distribution) is 
based on the well known relation 

1 py- T iirSix) (38) 



x ± it x 
Then, the spectral function turns out to be 



y(\) = -- lim lmG(z)\ z=x+le (39) 

7T e^O 
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In the case of the Gaussian Green's function (j37|), the discontinuities in the 
Green's function come from the cut of the square root, leading to the Wigner's 
semicircle law i'(A) = 57 V4 — A 2 for the eigenvalue distribution of random 
hermitian matrices. 

The Blue's function is defined as a functional inverse of the Green's func- 
tion 

B[G(z)] = z (40) 

For random Gaussian ensemble, the functional invepse of (^) is simply B(z) — 
z + 1/z. Blue's functions were introduced by ZeeEJ as the result of diagram- 
matic reinternretation of the seminal addition formalism (R-transformation) 
by VoiculescuB. 

3.1 Addition 

Let us consider the problem of calculating Green's function for the sum of two 
independent ensembles Mi and M2, i-e. 

CW = i J [dM 1 }ldM 2 ]P(M 1 )P{M2)tr - —J—^ (41) 

1 V 1 



N \ z - Mi - M% 

The concept of addition law for hermitian ensembles relies on additive trans- 
formation, which linearizes the convolution of non-commutative matrices (^) , 
alike to the logarithm of the Fourier transformation for the convolution of arbi- 
trary functions. This method is an important shortcut to obtain the equations 
for the Green's functions for a sum of matrices, starting from the knowledge 
of the Green's functions of individual ensembles_of matrices. 
The addition law for Blue's functions readsEl 

B x+2 {z)=B x {z) + B 2 {z)-~. (42) 

z 

Using the definition of the Blue's function we could rewrite the last equation 
in the "operational" form 

z = B X (G) + B 2 (G) - G- 1 (43) 

with G1+2 = G. The algorithm of addition is now surprisingly simplcS: Know- 
ing G\ and G2, we find ( |40| ) B\ and B2. Then we read out from j43| ) the final 
equation for the resolvent for the sum. Note that the method treats on equal 
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footing the Gaussian and non-Gaussian ensembles, provided that the ensem- 
bles are sufficiently independent (free). It is also applicable for chiral random 
matrices. For generalization of this algorithm-jar, the case arbitrary complex 
random matrices we refer to original papers EjlSeZI. 



3.2 Multiplication 

Blue's functions provide also an important shortcut to obtain the equation for 
the Green's (Blue's) function for a product of matrices, starting from knowl- 
edge of Green's (Blue's) functions of individual ensembles of matrices, i.e. to 
find 

G{z) = \- T (tr -j — — \ (44) 

w N \ z- Mx ■ M 2 / 

provided that Aii and M.i are free and (trA^i), (trA^) 7^ 0. Introducing the 
notation for continued fraction 

(45) 



■'M 



->M\ 



where self-energy S(z) is defined in a standard way 1/G(z) = z — S(z) multi- 
plication law for Blue's functions readH 

„ ( w \ w „ ( w \ „ ( w \ . 

b A^=—jA^))- b A^) (46) 

The usage of Blue's functions provides the factorization mechanism for aver- 
aging procedure in (Q) and is equivalent to S transformation of Voiculescu 
viaEj 

S(w) = / (47) 
where S(w) = x(w)(l + with \{zG{z) — 1) = l/z. 



3.3 Differentiation 

As demonstrated by Zee, differentiation of the Blue's function leads to deter- 
mination of the endpoints of the eigenvalue distribution for hermitian random 
matrices. One could easily visualize this point recalling the form of the Green's 
function for Gaussian distribution (|37]). When approaching the endpoints 
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z = ±2, the imaginary part of Green's function vanishes like y/{2 — A)(2 + A). 
Therefore, the endpoints fulfill the equation G'(z)\ z =±2 — °°j which may be 
used as the defining equation for unknown endpoints in case of more compli- 
cated ensembles. Since Blue's function is the functional inverse of the Green's 
function, the location of the endpoints could be inferred from a simpler rela- 
tion 

3.4 Integration 

Let us consider the partition function depending on a parameter z: 

Z N = (det(> - M)) (49) 

For finite N the partition function is a polynomial in z and has N zeroes 
(Yang- Lee zeroes): 

Zn = (z - zi)(z - z 2 ) ■ ■ ■ (z - z N ) (50) 

Taking the logarithm and approximating the density of zeroes by a continuous 
distribution we get 

log (det(z - M)) = J2 lo g( z - z i) = j P( z ') lo g( z ~ Z 'W (51) 

i 

We see that the density of Yang-Lee zeroes can be reconstructed from the 
discontinuities of the unquenched Green's function: 

d z log (det(j8 - M)) = (tv ■ det(z - M)\ (52) 

using e.g. Gauss law @. For z close to infinity the unquenched and quenched 
( |35| ) Green's functions coincide^. Moreover the unquenched resolvent is non- 
singular configuration by configuration and hence is holomorphic. Therefore 
it is exactly the functional inverse of the ordinary hermitian Blue's function. 
In general we may have several branches of Green's functions, which are de- 
termined by different saddle points contributing to (^). The discontinuity 
(cusp) and hence the location of Yang-Lee zeroes is determined by the possible 
contribution of two of them. From a \/N expansion 

logZ N =NF + F 1 + ^F 2 + ... (53) 



e In leading order in 1/N. 
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we see that two saddle points may contribute if 



ReF sp / = ReF° pJI (54) 

and F is determined by (|5^ ) 

F = J dz'G(z') + const (55) 

or equivalently 

E = zG- JdG B{G) + const (56) 

after integrating by parts. .Note that we have used the fact that z(G) — B(G) 
is just the Blue's function El of G, The integration using Blue's function is in 
most cases much simpler than the corresponding direct integration of G{z). 



3. 5 Mapping 

In case of non-hermitian random matrices, eigenvalues are complex. The 
average density of eigenvalues follows from the electrostatic analogy to two- 
dimensional Gauss law 

v{z,z) = -dsG{z,z) (57) 

Contrary to hermitian case, when eigenvalues condense on cuts, here they form 
two-dimensional nonholomorphic (non-analytic) "islands" . Surprisingly, using 
only the analytical structure of the Blue's functions, we can find the boundaries 
of these "islands" , while staying wholly within the holomorphic regionCj. 

Let us consider an example. First, we consider the case where a Gaussian 
random and hermitian matrix H is added to an arbitrary hermitian matrix M, 
then the case where a Gaussian random and anti-hermitian matrix iH is added 
to the same matrix M. In first case, eigenvalues are on the cuts (they are real), 
in the second case, form an island. In both cases the regions free of eigenvalues 
are holomorphic (here fulfill the Laplace equation). This implies that there 
exists a conformal transformation relating these two domains of analyticity. 
For the case considered here, Blue's functions for both cases are related via 

B iH+M (u) = B H+M {u) - 2u (58) 

Substituting u — > Gh+m(z) we can rewrite ( p8| ) as 

B iH +M [G h +m (z)]=z- 2G„+m {z). (59) 
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Let w be a point in the complex plane for which Ghj+m{w) = Gh+m{z). 
Then, using the definition of the Blue's function, we get 

w = z-2G h +m(z). (60) 

Equation ( |60| ) provides a conformal transformation mapping the holomorphic 
domain of the ensemble H + M (i.e. the complex plane z minus cuts) onto 
the holomorphic domains of the ensemble iH + M, i.e. the complex plane w 
minus the "islands" , defining in this way the support of the eigenvalues. Such 
construction could be generalized for other non-gaussian ensembles, due to the 
fact that addition of Blue's functions forms an abelian group. 
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